Quantum theory for spatial motion of polaritons in inhomogeneous fields 
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Polaritons are the collective excitations of many atoms dressed by resonant photons, which can 
be used to explain the slow light propagation with the mechanism of electromagnetically induced 
transparency. As quasi-particles, these collective excitations possess the typical feature of the matter 
particles, which can be reflected and deflected by the inhomogeneous medium in its spatial motion 
with some velocity. In this paper we develop a quantum theory to systematically describe the spatial 
motion of polaritons in inhomogeneous magnetic and optical fields. This theoretical approach treats 
these quasi-particles through an effective Schrodinger equation with anisotropic depression that the 
longitudinal motion is like a ultra-relativistic motion of a "slow light velocity" while the transverse 
motion is of non-relativity with certain effective mass. We find that, after passing through the EIT 
medium, the light ray bends due to the spatial-dependent profile of external field. This phenomenon 
explicitly demonstrates the exotic corpuscular and anisotropic property of polaritons. 

PACS numbers: 03.65.-w, 42.50.Ct, 42.50.Gy 



I. INTRODUCTION 

Quasi-particles are excitations of the matter. Accord- 
ing to modern many body theory, Elementary parti- 
cles and quasi-particles are basic construction of mat- 
ters. The latter are crucial for understanding many phe- 
nomena in condensed matter physics. Actually, quasi- 
particles can be regarded as collective excitations of many 
elementary particles, as well as the mixtures of different 
elementary excitations, whose behavior are similar to the 
matter particles 

In atomic physics and quantum optics, some exotic 
phenomena can be explained by the concept of quasi- 
particles. For example, slow light phenomenon [i 3| 
in electromagnetically induced transparency (EIT) 0, [sj 
can be explained in terms of quasi-particles - polaritons 
S S i, li- EIT happens when a weak signal light field 
and a stronger control field are coupled to an ensemble 
of atoms with a A energy level configuration. Under the 
two-photon resonance, due to the destructive interfer- 
ence between two interaction paths, the initially opaque 
resonant medium becomes transparent with respect to 
the probe field, and the group velocity of light is slowed 
down. Light then is stopped in the EIT medium because 
only the dark state polariton is excited. The dark state 
polariton is a bosonic like collective excitation, which is 
a mixture of a signal light field and an atomic spin wave 

Most recently, the light deflection was observed for the 
EIT atomic medium in an external field with spatially 
inhomogeneous distribution [11, , ,12.] . In the experiment 
of Ref. [ll[ , it is found that the light ray bends when a 
magnetic field with small gradient vertical to the prop- 
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agation direction is applied to a cell with A-type rubid- 
ium gas. This experiment was interpreted as the Stern- 
Gerlach experiment of the dark polariton, thus the effec- 
tive magnetic moment of the dark state polariton is ob- 
served for the first time. It demonstrates that the dark 
state polariton indeed behaves as a matter particle with 
mass, momentum and magnetic moment et al, which can 
be reflected, refracted, and even deflected by a gradi- 
ent force. Therefore, quasi-particles show their particle 
nature with definite momentum and effective mass. Dif- 
ferent from that in Ref. [11], the experiment of Re. (l^ 
shows that a light can also been deflected by an optical 
driven Rb atomic vapor when the profile of the driving 
field is inhomogeneous. In this situation the angle of de- 
viation is an order of magnitude larger than that in Ref. 
[Tl| . The observed phenomenon about light defiection in 
such EIT media has been explained correctly according to 
the semi-classical theory [l^ without using the concept 
of dark-state polariton, which needs the quantization of 
light fields. 

Like matter particles, quasi-particles possess the wave- 
particle duality, that is, quasi-particles sometimes appear 
to behave as particles and sometimes appear to behave 
as waves. Here, we are interested in the particle aspect 
of the dark polariton, which is an atomic collective ex- 
citation dressed by the quantized probe light. The main 
purpose of this paper is to systematically develop a quan- 
tum theory describing the spatial motion of polaritons in 
inhomogeneous magnetic and optical fields. We begin 
our investigation with the propagation of quasi-particles 
in the limits of atomic linear response, where the atomic 
equations are treated perturbatively. With an effective 
potential induced by the steady atomic response in the 
external spatial-dependent field, the dynamics of spatial 
motion of the quasi-particles is governed by the effective 
Schrodinger equation. The spatial motion of the quasi- 
particle is of anisotropic depression - the longitudinal 
motion is like a ultra-relativistic motion of a "slow light" 



while the transverse motion is of non-relativity with cer- 
tain effective mass. 

This paper is organized as follows: In sec. |TT1 we 
present the theoretical model for a A-type atomic en- 
semble in the presence of inhomogeneous external fields, 
and derive the system of equations governing the spatial 
motion of the signal field in the atomic linear response 
with respect to the probe field. In Sec. Illli the per- 
turbation theory is applied to obtain the atomic motion 
equation which is related to the linear response to the 
signal field. In Sec. IIVI and Sec. |Vl the crucial idea of 
the EIT - the dark-state polariton is introduced as an 
dressed fields to describe the spatial motion of collective 
excitation. Afterward, the dynamics of the quasi-particle 
- dark polariton is discussed in the presence of an inhomo- 
geneous magnetic field with a spatial distribution along 
the transverse direction. In Sec. IVIl the spatial motion 
of the signal light in an inhomogeneous coupling field is 
investigated. Then we make our conclusion in Sec. IVIII 



II. THEORETICAL MODEL FOR A-TYPE 
ATOMIC ENSEMBLE IN EXTERNAL FIELDS 

We consider an ensemble of N identical and noninter- 
acting atoms, which is confined in a cell ABCD as shown 
in Fig. [T|d. Each of the atoms is modeled by a A-shaped 
energy level configuration with internal states jg), \s) and 
|e). The transitions from the two lower states \g) and |s) 
to the excited state |e) are coupled by two optical fields, a 
weaker probe field and a stronger control field, as shown 
in the top panel of Fig. [1^. The atomic transition from 
\g) to |s) is forbidden by the electronic dipole coupling. 
The probe field carries frequency f and the wave num- 
ber k. It is a quantized electromagnetic field with (t+ 
polarization. Under the rotating wave approximations, 
its negative frequency part of the electric field (r, t), 
couples the ground state \g) to the excited state |e) at res- 
onance, in the absence of the magnetic field B (r) . The 
control field has carrier frequency Vc = lo,, — lOs and wave 
number kc- It is a classical field with cr~-polarized, and 
couples to the upper state |e) and the metastable state 
\s) with Rabi frequency ^l (r). After the magnetic field 
is applied along the z-direction, the internal energies of 
the corresponding states are shifted from their origins by 
magnitudes fXiB with 



pij = mpgpfiB,i G {g,s,e}. 



(1) 



Here, fis is the Bohr magneton, gp is the Lande g- factor 
of the internal state i, and m^^ is the magnetic quantum 
number. 

As shown in Fig. [TJd, the probe field and the control 
field propagate parallel in the z— direction with wave 
number k and kc respectively. The Hamiltonian of this 
typical EIT system is given by iJ = + if + . 
Let us use ct^^ (t) = to denote the internal state 

operator of the j-th atom between states and {ly). We 




n —WW 



Signal field 




Control field 



FIG. 1: (Color on line) (a) Level scheme of atoms interact- 
ing with (T"*"-polarized probe and (T~-polarized probe strong 
fields. Q denotes the Rabi frequency of the undepleted classi- 
cal control field, (b) Configuration of the optical beams and 
the magnetic field inside the atomic medium. 



introduce the collective atomic operator Q 



C^M^ (r, = ]^ X! 



(2) 



which is averaged over a small but macroscopic volume 
containing many atoms Nr — [N/V) dV ^ 1 around 
position r. Here N is the number of atoms in an inter- 
action volume V . Then the Hamiltonian of the atomic 
part reads 



N ^ 



(3) 



where we have neglected the kinetic term of atoms, and 
Uj are the corresponding energy level spacing of the inter- 
nal atomic level. iJ^^^ is the free Hamiltonian of the radi- 
ation field. Using the electric-dipole approximation and 
the rotating-wave approximation, the interaction with 
electromagnetic field reads @, 0, H, |^ 



N 
V 



h.c. 



(4) 

Here, is the negative frequency of the probe field; 
f2 (r) is the Rabi frequency of the control field, which 
usually depends on the spatial coordinate through the 
spatial profile of driving field; d^g is the dipole matrix 
element between the states \g) and |e). 

For convenience, we describe the electric field as 



E+ {r,t) 



E{r, t) e*(fe^-'^*) 



(5) 
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in the following discussion. Here exp [i {kz — vt)] is the 
carrier wave with frequency v and wave number k propa- 
gating in the z direction, and E (r, t) is the slow varying 
envelope, meaning that its spatiotemporal variation is 
much slower than the carrier wave length and frequency. 
Further we introduce the slowly varying variables for the 
atomic transition operators 



CTeg (r, t) = CTeg (r, t) e 
CTes {r, t) = CTes (r, t) e 



-ikz 



(6a) 
(6b) 



In the rotating reference frame, the dynamics of this sys- 
tem is described by the interaction Hamiltonian 



N f 

H, = -- d\^^l,B<7,J 



(7) 



-I- {gaegE + Qaes + h.c.)], 
where the atom-field coupling constant g is defined as 



2eoV 



(8) 



Before we study the EIT features of this system in de- 
tail, let us first stand in the point of view of light to 
investigate the propagation effects of pulses in an atomic 
medium. It is well-known that, when atoms are sub- 
jected to an electric field, the applied field displaces the 
positive charges and the negative charges in atoms from 
their usual positions. This small movement that positive 
charges in one direction and negative ones in the other 
will result in collective induced electric-dipole moments. 
All dipole moments in the dielectric material generate the 
polarization collectively, which is defined as the collective 
dipole moment per unit volume 



P 



N 

r 



J{kz- 



h.c. 



(9) 



The collective dipole moment in Eq.(l9]) is caused by the 
atomic response to an optical electric field in a dielectric 
material. In turn, every dipole with a nonvanishing sec- 
ond derivative in time radiates an electromagnetic wave, 
that is, the dielectric response P of the medium acts as 
an effective source to produce the electromagnetic field. 

The Heisenberg equation for the slowly varying field 
operator E (r, t) results in a paraxial wave equation in 
classical optics [1, [3| 



. d . d c 2 

'dl^'^d-z^Yk^^ 



E = ~^g*Na^ 



ge- 



(10) 



Here, c = l/^eoMo is the velocity of light in vacuum and 
the transverse Laplacian is defined as 



(11) 



in the rectangular coordinates. When we neglect the x- 
and y-dependence of E, that is, confine the problem in 



one dimension, Eq. ( I10[) immediately reduces to the 
usual propagation equation 



d_ 

di 



d_ 

dz 



E 



~9*Nag 



(12) 



given in Ref. [1, 0, H, H, [T3|i which only describe light 
propagation in z direction. To consider the problem in 
three spatial dimensions, one can use the paraxial wave 
equation ( [10] ) to investigate the dynamics of the input 
pulse in a resonant atomic medium. 

In this paper we will focus on the case that the linear 
optical response theory works well, which can sufficiently 
refiect the main physical features of the spatial motion of 
the input pulse with slow group velocity. The lowest or- 
der contribution to the polarization is the linear response 
of atoms defined as 



p(i) 



N 



(I) i{kz-u 
eCTge e 



(13) 



Then, the paraxial wave equation becomes 

where the definition of a'^ge will be given in the next sec- 
tion. 



III. 



PERTURBATION APPROACH 



We now study the evolution of the atomic ensemble 
under the influence of the applied optical fields. The 
dynamics of this atomic ensemble is described by the 
Heisenberg equations 

c^eg = [i (Mff - ^^e) B] (T^g - in*asg (15a) 

+ ig* ((Tee - CTgg) E^ , 

(7es = [i (Ms - A^e) B] a^s - ig*agsE+ (15b) 

-t- ifl* ((Tee - CTss) , 
&sg = [i (Ms - Ms) B] (Tsg - ifl<Jeg + ig*<JseE^. (15c) 

Since EIT is primarily concerned with the nonlinear mod- 
ification of the optical properties of the probe field, thus 
the low density approximation is valid. In this approxi- 
mation, the intensity of the quantum field is much weaker 
than that of the coupling field fl, and the number of pho- 
tons contained in the signal pulse is much less than the 
number of atoms in the sample. 

In the low density approximation, the perturbation ap- 
proach can be applied to the atomic part, which is intro- 
duced in terms of perturbation expansion d, 0, H 111 



,(0) 



(16) 



where 



{e,s,g} and A is a continuously varying 
parameter ranging from zero to unity. Here a^^^ is of the 
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zeroth order in gE, is of the first order in gE and 
so on. We now substitute Eq. into Eq. and 

retain only terms up to the first order in the signal field 
amplitude. We thereby obtain the system of equations 
in the zeroth order 



^(0) _ J* (0) 



eg 



(0) 



.(0) 



where the parameters 



di = i - fig) B - 71, 
d2 ^iifJ,s- t^g)B ~ 72, 
dz ^i{fi^- /is) B - 73. 



(17a) 
(17b) 
(17c) 



(18a) 
(18b) 
(18c) 



and we have phenomenologically introduced the energy- 
level decay rates 7i( i £ {1, 2, 3}). 

We assume that all population of atoms are initially 
prepared in the ground state \g) in the absence of elec- 
tromagnetic fields, and the depletion of the ground state 
is not significant for any time t > due to the quantum 
interference effect, therefore 



= 1 



(19) 



while others vanish ^21- Then, the first order atomic tran- 
sition operator crjj\ which are related to the atomic lin- 
ear response to the probe field, satisfy the following equa- 
tions 



eg 



dlag') -tn*ag'> -tg*E\ (20a) 
cri'J ^d;cTil^+^n*[ail^-ag^), (20b) 

(20c) 



"sg — "'2"sg '"^^eg 



In order to get the equations of motion for polaritons, we 
rewrite Eq. (|^D|) as 



{dt~di)-{dt-d2) + n 

idt-d2)a<^'K 



gs ' 



(21a) 
(21b) 



Equations ([T^ and (PT|) constitute a self-consistent 
system of equations, which indicates that the polariza- 
tion field age"* can serve as a source to generate the elec- 
tric fields, whereas the propagating light in turn drives 
the atomic media via the dipole interaction. They are 
the starting point of our investigation in the following 
several sections, where we study the phenomena of light 
deflection that occur as a consequence of the interaction 
between the A-type atomic ensemble and an external field 
with a spatial distribution. 



IV. SPATIAL MOTION OF QUASI-PARTICLE 
IN A HARMONIC MAGNETIC FIELD 



It is well known that the EIT system has two remark- 
able properties: 1) the opaque absorption medium be- 
comes transparent with respect to the probe light at cer- 
tain frequencies. It happens because the absorption on 
both transitions is suppressed by the destructive interfer- 
ence between the excitation pathways to the upper level. 
Thus a transparency window is rendered over a narrow 
spectral range within the absorption line. 2) The group 
velocity of the incoming pulse has been largely reduced 
within the transparency window. Physically, the slow 
light in a EIT system is interpreted by the formation of 
so called dark-state polariton (DSP). A "dark-state po- 
lariton" is a bosonic-like collective excitation of a signal 
light field and an atomic spin wave [1, 0, H, , whose rel- 
ative amplitude is determined by the control laser field. 
In this section, we study the dynamic of the DSP in the 
presence of a harmonic field with a spatially inhomoge- 
neous distribution in the transverse direction, where the 
control field is assumed to be independent of position and 
time. 

When the light pulse enters a medium, photons in- 
teract with atoms of the medium. They then combine 
together to form a type of excitations known as polari- 
tons, which are one kind of quasi-particles. In an EIT 
system, two types of polaritons are introduced - the dark 
polariton and the bright polariton, which are described 
respectively by the dark polariton field operator ^E* and 
bright polariton field operator $ 



* (r, t) ^ Ecos0 
<^{r,t) ^ Esm0^ 



iVCT<,V sin ( 



(22a) 
(22b) 



They are atomic collective excitation (quasi-spin wave) 
dressed by the quantized probe light with the inverse 
relations 



iJ = * cos ( 



<i>sin( 



1 



- ^= ($ COS 6* - sin 6*) 



(23a) 
(23b) 



The dark polariton field operator ^E* and bright polariton 
field operator $ have bosonic commutation relations in 
the limit of few photons and many atoms. And the ac- 
tion of ^'t on the vacuum creates the dark states, which 
contain no component of the excited state |e). By assum- 
ing that the Rabi frequency is real, the mixing angle of 
the signal field and the collective atomic polarization is 
given by 



tan 9 = 



(24) 



where the Rabi frequency il is related to the control laser 
power P through fl^ — 2 \des\^ P/ {ceoS). In the nearly 
two-photon resonant condition, the only excitations are 



5 



dark polaritons, which generate an eigen-state with van- 
ishing eigenvalue of the interaction Hamiltonian. It is 
found from Eq. ([22|) and ( [24]) that, by reducing the am- 
plitude of the control field, the contributions of light or 
atoms to the DSP can be changed, then the DSP varies 
from photons to atoms. Thus it is roughly seen that 
the mixing angle 9 determines whether or not the group 
velocity of the signal pulse propagating in the atomic 
medium can be decreased. 

In order to find how the mixing angle affects the group 
velocity of the input pulse, we derive the equations of the 
spatial motion for the dark polariton fields. In the above 
sections, we have achieved the dynamic motion equations 
of atoms and light. In terms of the field operators for the 
dark and bright polaritons, Eq. ([2T|) and ( [T4|) can be 
rewritten as 



(25) 



and 



gVN {-^ cos 6* + $ sin 9) = 

- [{dt - di) i {dt - d2) + n] 
X (^cos^ - ^fsin^) 



{dt - d2) ($cos6'- *sin^ 



gVN 
in 



(26) 



For a very small magnetic fields, we have |/Ltg — /Ze| <C 71. 
Furthermore, we assume a sufficiently strong driving field 
such that fl^ 3> 7i72. In the adiabatic approximation, 
the excitation of the bright polariton field $ vanishes 
approximately. Then the dynamics of the dark polariton 
field ^' is governed by the Schrodinger-like equation 



with an effective potential 



(27) 



(28) 



induced by the steady atomic response in the external 
spatial-dependent field. Here, we have set 72 = 0; while 
the effective kinetic operator 



f = vgp, - ^ cos^ 



(29) 



represents an anisotropic depression, where the momen- 
tum along z direction is defined as Pz = —id^- The longi- 
tudinal term VgPz in Eq. (|29p describes a ultra-relativistic 
motion with a slow light velocity 

Vg^ccos'^9, (30) 

while the transverse part Py/(2m) in the effective kinetic 
term desribes a non-relativistic motion with an effective 
transverse mass 



m = — = — scc^ 
v„ c 



(31) 



The above effective Schrodinger equation governs the dy- 
namics of spatial motion of quasi-particles. 

Obviously, when no magnetic field is applied, due to 
the transverse Laplacian operator V|i commutating with 
dz, we can separate the z-component from x-y compo- 
nent. Neglecting the x- and y-dependence of 4*, Eq. P?)) 
describes a stable propagation along z- axis with group 
velocity Vg [1, H, H, 0, fi @l- Hence, the amplitude of the 
control field determines the group velocity of the input 
pulse in the atomic medium. Adiabatically rotating the 
angle from to tt/2, the polariton can be decelerated to a 
full stop. On the other hand by increasing the strength of 
the coupling field, that is, reversing the rotation of 9 adi- 
abatically, it leads to a re-acceleration of the dark-state 
polariton associated with a change of character from col- 
lective spin-like waves to electromagnetic photons. 

Now we consider the three-dimension problem. By 
defining Pj = —idj{ j € {x, y, z}), Eq. lfST]) can be rewrit- 
ten as an effective Schrodinger equation 



with the effective Hamiltonian 



H, 



eff 



2m 



(32) 



(33) 



where /i' = (/is — fJ-g) sin^ 9. The magnitude of the effec- 
tive transverse mass is totally determined by the mixing 
angle 9 of the signal field and the collective atomic polar- 
ization. When the amplitude of the control field is small, 
spin waves have large contributions to the DSP, therefore 
the effective transverse mass is large; when the Rabi fre- 
quency is large, photons give large contributions to the 
DSP, therefore the effective transverse mass is small. The 
effective Schrodinger equation (|32p is the starting point 
for investigating the spatial motion of the dark polariton. 
It shows that, due to the inhomogeneity of the magnetic 
field, the motion of the dark polariton will be scattered 
by an effective potential with value 11' B (r). 

Now we assume that the magnetic field in z direction 
has a spatial distribution in the transverse direction with 
the expression 



B (r) - Bo + B^x^ + Byv^ 



(34) 



where B^, By < 0. Then the effective Hamiltonian oper- 
ator becomes Hi = H„ia + VgPz — fJ-'Bg with 



2™ 



2to 



(35) 



In classical physics, H„ia corresponds to a two dimen- 
sional harmonic oscillator with mass m and angular fre- 
quency ujx = —2fi' Bx/m in x-direction and ujy = 
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yJ—2^'By/m in y-direction. For a given initial state 
^ (0), the evolution state 4" (t) of the system is a unitary 
transformation of the initial state ^ (0) with the time- 
evolution unitary operator U [t) = Uz (t) Uy {t) Ux (t) : 



Uz = cxp {-ivgPzt) 
Uy = exp 



Ux = exp 



2m 



a: 2 1. 

—a:; t 



(36a) 
(36b) 

(36c) 



Next we consider the evolution dynamics of a spa- 
tially well-localized wave packet, which is centered at 
(xQ,yQ^ zq) — (0,0,0) and has a vanishing mean veloc- 
ity in all directions. The spatially well-localized wave 
packet is assumed to be initially in a Gaussian form 



*(o)= n 



1/4 



(37) 



with width l/y^o^, ^ S z} in a;, y, z-direction, re- 

spectively, and «(; 7^ Xj where Xj = mujj. By getting rid 
of an irrelevant global phase factor, the wave function at 
time t > reads 



I TT 



1/4 



(38) 



E 

Tiin2 



The initial Gaussian wave packet evolves into a superpo- 

mtum num- 
Cocfficicnts 



sition of the product states (f)!^^ (j)n2 with quantum num- 



bers ni and n2 taking on the value 0,1,2- 
in Eq.dSll) read 



^2n — 



Here, (j)'n'^ is the eigenfunction of Hamiltonian P^/ (2m) 
maj^x^/2 with the corresponding eigenvalue En 

{n + l/2)ujx 



2"n! 



(A, 



1/4 



Aj - 

A, -I- an 



(39) 



2"r7,! 



A. 



(40) 



where {x) are Hermite polynomials. The wave func- 
tion ^i^^ has the similar expression as Eq. ipD)) with x 
replaced by y. 

Since the wave-function at time t is an even function 
of variables x and y, which can be seen from Eq. (j38p 
to Eq. the trajectory of the center of the wave 

packet in the x - y plane does not change with time. But 
in z-direction, the dark polariton propagates with mean 
velocity Vg, and the center of the wave packet leaves its 
original place proportionally to the time t with (z) = Vgt. 




FIG. 2: (Color on line) (a) Schematic illustration about the 
time evolution for the center of the initial Gaussian state of 
the dark polariton along z direction, (b) the time evolution 
of the variance in a;-direction. The variance is in unit of 0.1. 



Although the wave packet keeps it shape along z direction 
with an unchanged variance 



(Az)^ 



1 

4a^ 



(41) 



the variances in the x- and y-direction oscillate with time, 
namely, 



(Ax) = A^xCos{2ujxt) + A+x, 
{Ayf = A^y cos {2ujyt) + A+y, 



(42) 



where 



A 



±s 



2'Km?uj1 ± at 
SagTrm^Wg 



-,s = x,y. 



(43) 



Fig. [5] (a) schematically illustrates the time evolution of 
the initial Gaussian packet. The wave packet distributed 
along the z-direction keeps its original shape. However, 
in x-direction, the shape of the Gaussian packet oscillates 
as time evolution, its width change is shown in Fig. [2](b). 

In plotting Fig. [2] (b) , we have taken the reasonable pa- 
rameters accessible in current expriments pT| : the width 
of the initial Gaussian 1/ax = 1 /4cm, the atomic density 
N/V = 10^^ per cube centimeter, the quantum num- 
ber of the ground state mPp = —2 with g^p = 1/2, the 
quantum number of the mestable state rrip = 0, and 
71 = 27r X 2.87MHz. 



V. THE DEFLECTION OF DARK 
POLARITONS IN A LINEAR MAGNETIC FIELD 

In a recent experiment [Tlj . a magnetic field with 
small transverse gradient is applied to a A-type atomic 
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medium. It was observed that the Ught beam is deflected 
after the signal light passing through the EIT gas cell. 
This observation first demonstrates that quasi-particles 
- dark-state polaritons have a non-zero magnetic mo- 
ment. This experimental observation can be interpreted 
straightforwardly according to the above quantum theory 
of spatial motion of polaritons in inhomogeneous fields. 
For simplicity, we assume the magnetic field 



B (r) =Bo + Bix. 



(44) 



has a linear gradient along the x-direction. We further 
allow the input pulse to vary only in one transverse di- 
mension, says, in the x-direction, which means that we 
neglect the y-dependence of the input pulse E. Then 
the two dimensional effective Hamiltonian for the dark 
polaritons reads 



Ho 



VgPz 



Im 



fiCx, 



where the parameters 



C = 2Bi sin^ 



bo 



2Bo sin^ 0, 



(45) 

(46a) 
(46b) 
(46c) 



can be controlled by the mixing angle 9. For an initial 
dark polariton field with Gaussian spatial distribution 
Oix = otz = as given in Eq. (|37p . the time evolution 
of the polariton field reads Uu (t) = exp {—iH2t). By the 
Wei-Norman algebraic method (see the appendix A) [31 , 
the unitary operator Uu (t) can be factorizcd as Uu (t) = 
U2 (t) Ui (t) with 



U2 it) = e 
Ui{t) 



_p2 if2_MC.p 



(47a) 
(47b) 



A straight forward calculation shows that, the initial 
Gaussian packet evolves into 



I/tt 



„Vt(fco-V^ 



(48) 



e*^^*^ exp 



\ Ira j \ m / 

264 + 2t2/m2 



At time t, the center of the input pulse moves to Vgt along 
the z-direction and ^C, / {2rri) along the x-direction. 
When a dark polariton is excitated by the interaction 
between light and atoms, the dark polariton will achieve 
a velocity along the x-direction as 



mvo 



(49) 



after it pass through the gas cell with length L. Therefore 
the deflection angle reads 



L fi 

Vgk 



Bi sin^ 9. 



(50) 




FIG. 3: (Color on /mejSchematic illustration about the spa- 
tial resolution. 



In real experiment, the dephasing time is nonzero due 
to the collision between atoms, which leads to the ab- 
sorption of the energy of the probe beam by the atomic 
medium. 

The above results mean that the deflection angle of 
the output pulse depends on the mixing angle 9 between 
the signal field and the collective atomic polarization, the 
wave number k of the input pulse, and the gradient Bi of 
the inhomogeneous magnetic field. One can find that the 
magnetic moment of the dark polariton has an effective 
value 



Mpo/ = Aisin^ 9. (51) 
-2 and m, = 0, we find the effective 



By taking rUg — - 
magnetic moment 



2g^^\iB sin^ ( 



(52) 



which is exactly the theoretical result given in Ref . [Uj . 

Next we consider the spatial resolution, which in optics 
refiects the ability of this optical system to form separate 
and distinct images of two objects. The spatial resolution 
is defined here as the mean signal divided by its standard 
deviation 



R 



Ax 



62 



2m264 + 2t2 ' 



(53) 



where the mean position (x) in the transverse direction 
and the standard deviation Ax are given by 



+00 



(t) x^ {t) dxdz = 



2m' 



Ax 



'x2) - (x)^ 



m264 + ^2 



(54a) 



(54b) 



It can be seen that the spatial resolution increases as the 
interaction time between light and atoms increases. 

Actually, the phenomenon of light deflection in such an 
inhomogeneous magnetic field can also be described with- 
out using the concept of quasiparticles - dark polaritons. 
Here, we show how to calculate the deflection angle a in 
Eq. ([50]) according to the semiclassical theory. We begin 
the semiclassical approaches by considering the evolution 
of the system from Eq. (fT4|) and ((20)) . First, the atomic 
linear response to the signal field has been explicitly re- 
fiected in Eq. (I20p . Under the adiabatic approximation 
that the evolution of the atomic system is much faster 
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than the temporal change of the radiation field, we can 
obtain the steady-state solution for the atomic transition 

(J^e by setting all time derivatives to zero in Eq. ([20]) . 
namely 



[« il^g - IJ.s)B + 72] g 



did2 



\n\ 



E 



(55) 



Here, the undepleted control-field approximation is used 
and 72 ~ is assumed. This approach based on the 
atomic linear response results in an effective potential for 
the motion of signal slow-varying amplitude due to the 
spatial distribution of the magnetic field. The spatial 
motion of the slow varying amplitude is described by the 
following equation: 



-fiBE 



(56) 



which describes a shape-preserving propagation in z di- 
rection with velocity c. 

For an initial Gaussian wave packet of i? in a; — z plane, 
after passing through the gas cell, the wave center shifts 
from {x,z) = (0,0) to the well-defined position 



(a;,z) 



2kc 



(57) 



And it can be obviously found that Eq. ([50|1 is the de- 
fiection angle in this approach. 



VI. DEFLECTION OF LIGHT IN 
INHOMOGENEOUS COUPLING FIELD 

In this section, we turn our discussion to the deflection 
of slow light in the atomic medium driven by a optical 
field with an inhomogencous profile, while the magnetic 
field is uniform. This phenomenon was experimentally 
observed in Ref. [13] where the cell filled with EIT-based 
atomic gas is referred as an ultra-dispersive optical prism 
with an angular dispersion. We note that the probe light 
is relatively strong in comparison with the control light in 
this experiment, thus the susceptibility obtained from the 
linear response theory can not work well to explain the 
experiment phenomenon. In this paper, we do not plan 
to explain the experiment data in Ref. [12l | in the strong 
coupling limit. Our main purpose is to predict a new 
quantum coherent phenomenon for the light deflection 
by the atomic media when the experiment is carried out 
for the weak probe field. 

We assume the strong driving field has a Gaussian pro- 
file 



O = f2o exp 



2(72 



(58) 



in the transverse direction. Here, we confine the problem 
to two dimensional space, the x-z plane. Then the trans- 
verse Laplacian operator reduces into a one dimensional 
operator — d'^/dx^. 

By invoking the steady-state conditions, it is found 
that the polarization field CTge^ , which serves as a source 
for the electric fields in Eq. is proportional to the 

slow varying amplitude E given in Eq. (|55p. Under a 
strong, undepleted driving field approximation, the cou- 
pling between atoms and light induces a spatial depen- 
dent potential into the propagation equation. The spatial 
shape of this potential induced by CTge'' is completely de- 
termined by the profile of the Rabi frequency which 
can be seen in the first line of Eq. Thus, when the 

signal pulse parallel to the control beam travels across 
the atomic cell, it will be scattered by the effective po- 
tential. However, as the width of probe beam is less than 
that of the control beam, the trajectory of the probe light 
may bend when it is adjusted to the left side or to the 
right side of the control beam profile; hence the probe 
and control beams are no longer parallel after they go 
through the gas cell. 

In order to investigate this phenomenon, we assume 
the probe beam is in a Gaussian state 



E{0,x,z) 



1 



: exp 



{x — a) 



z 

2&2 



(59) 



before it enters the gas cell, where b{< a) is the width of 
the probe field and a is the initial location of the wave 
packet center of the probe fleld along x-direction. The 
sign of a indicates the incident position comparatively to 
the left or right hand side of the control beam's center 
xo = 0, and the magnitude \a\ denotes the distance from 
the control beam's center. In order to investigate the 
evolution of this initial state, we expand jfip at the 
position a and retain the linear term proportional to x—a. 

With the above considerations, the paraxial equation 
in Eq. ([5S|) becomes 



iE + icd.E + ^dlE = (770 + m^) E 



where 



-"o-Ml-|^) M^^Aexpfe 



-2aA 



\9\'N 



^0 exp 



(60) 

(61a) 
(61b) 



and A = (^g — Hg) B. By making use of the Wei- Norman 
algebraic method [l3|, it is shown that, after passing 
through the Rb gas cell, the center position (x, z) — (a, 0) 
of the probe field at time t = is shifted to 



X = a -\ 
z^L. 



a^kc ' 



(62a) 
(62b) 



9 



z 





—1 ; 


: \ 
, \ 




• 
• 

• 

■ 


1 : 

1 ; 
> 

r : 
^^^^ 


1 »^ 
1 t 

■ t 

■ 1 
I 1 
• 1 

L 1 





FIG. 4: (Color on Zine^Schematic illustration about the ray 
deflection of the probe light in the presence of inhomogenous 
coupling light. The yellow line is the spatial distribution of 
the control light. 



If we track the motion of the center of the probe beam, 
a mirage effect occurs. The sign of A and the incident 
position a of the signal light determine whether the tra- 
jectory of probe beam bend. When the magnetic field is 
absent A = or the center of the probe field is coUinear 
to that of the control field a = xq = 0, the trajectory 
of the signal light is a straight line. We assign the posi- 
tive sign for a as the probe beam is shifted to the right 
with respect to the center of the control light, and denote 
a < as the signal beam is shifted to the left. When the 
probe beam is shifted to the right, in the case of A < 0, 
the signal light feels a "repulsion potential" due to the 
coefficient 771 > 0, thus the trajectory bends to the left; 
at the condition A > 0, the signal light undergoes an "at- 
tractive potential" in the atomic medium due to 771 < 0, 
thus the trajectory bends to the right. When a < 0, it 
can be found from Eq. (pT|) that due to the coefficient 
of the linear potential larger than zero, i.e. A > 0, the 
probe beam experiences a "repulsion potential" within 
the EIT medium, and its center is shifted to the left. As 
rji is smaller than zero, i.e. A < 0, the probe beam suffers 
an "attractive potential" during its passing through the 
EIT medium, hence its center is shifted to the right. The 
corresponding schematic diagram is given in Fig|4l where 
yellow solid line is the spatial distribution of the control 
light, the red dash lines give the deflection at A < 0, the 
blue dotted lines describe the light trajectory at A > 
and the black solid lines depict the light ray at A = 0. 
The same results about the deflection of light ray have 
been discovered by us using the semiclassical theory [l^ . 

From the point of particle nature, the force acting on 
the particle is completely determined by the value and 
sign of r]i . Thus for a particle passing through the point 
a 7^ 0, when A = 0, the particle does not feel any force, so 
it travels across straightly, so does the particle at point 
a = 0. For a particle traversing through the position 
a (> 0), when A < 0, this particle is subject to a negative 
force, which moves the particle to the left with respect to 
it original place; however, when A > 0, it experiences a 
positive force, which makes the particle move to the right. 
Also for a particle going through the place at a(< 0), 
when A < 0, the particle moves to the right because of 



the action of a positive force; when A > 0, it goes to the 
left due to the action of a negative force. 

We have to point out that the magnetic field is not 
necessary for the occurrence of the above-described phe- 
nomenon. For the model containing a A-type atomic en- 
semble interacting with one control beam and one probe 
beam, similar phenomenon can also be found as long as 
the two photon detuning 

A = Ap - Ae (63) 

varies, where Ap = tOeg — v is the detuning between the 
atomic transition from |e) to \g) and the probe beam, 
Ac = ojes — is the detuning between the atomic tran- 
sition from |e) to |s) and the control beam. In order to 
clarify the dependence of light deflection on two photon 
detuning A, we begin our description with the Hamil- 
tonian in the interaction picture. In the rotating frame 
with respect to 

Hq = UJe<7ee + (^^ + A^) (Jss + [uJg + Ap) CTgg , 

the interaction Hamiltonian reads 
N f 

Hj = -— / d^r[Apagg+Acass + {gcregE + flcFes + h.c.)], 

in the absence of magnetic field. The first order atomic 
transition operators have the similar form as Eq. ((20l) by 
replacing /x^B, ^gB and by Ac, Ap and zero re- 
spectively. Then the atomic transition operator ag]} = 
gAE/ \ ^l\^ , which induces a potential dependent on the 
two photon detuning A. 

In a real experiment, the dephasing rate of the forbid- 
den |e)-|s) transition is nonzero due to atomic collisions 
etc. Therefore an additional anti-Hermitian decay term 
will be introduced phenomenologically into the effective 
Hamiltonian 

Heff = cp^ + ^pI + v'q + ?][ {x ~ a) , (64) 

where r]'^ = —aj — ihj, j = 0, 1 are complex. Then it can 
be found that, after light passing through the Rb gas 
cell, the dephasing rate introduces two additional terms 
to Eq. dSli) 

where T — L/c is the time for the light travelling through 
the medium, and 

2aN\g\^ fa^\ ^ 

For an atomic medium with length L — 7.5cm and den- 
sity N/V = lO^^CTO"'^, the dephasing rate 72 « 10~''7i. 
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When a control beam with width cr = L/A and fre- 
quency = 5 X 10^'* iJz is coupled to the atomic en- 
semble with D,Q = 571, for the probe beam with width 
0.07mm incident at position a = cr/2, two additional 
terms b'^biT « 10"^ ^nd bic^T^/{b'^P) « lO''^ caused 
by dephasing are much smaller than that of the term 
aicT^ / {2k) « 5 X 10~^ induced by the frequency detun- 
ing A, which means Eq. (|62b.) dominates the distance in 
the transverse direction. 

Next we investigate how the trajectory of the probe 
beam behaves when an effective potential include the 
quadratic term of x and y in the transverse direction. 
This induced potential is obtained when we expand 



I f2 1 ^ = exp 



2 9 
X y 



(66) 



around the center a.^ and ay of the incident beam with 



the profile shape 



E{Q) = 



exp 



z^ + (a:-Qx) + [y - ayY 



(67) 



By retaining the quadratic term of a; — a^; and y — ay. 
The paraxial motion equation becomes 

idtE + icd,E+-^{dl + dl)E= [V (x) -I- V {y)] E (68) 



where 

y (x) = [Cxo + Cxi (X - ax) + Cx2 (X - ay,y 
The coefficients for x — {^^ y\ ^-re 

-f^o"^ exp 

2\g\^N 



E. 



Cxi 



Cx2 



^5 

tI + 2a 



-flyAilQ ^ exp 



2 

Xr.-2 



l^o^lsl'^A^Aexp 



(69a) 
(69b) 
(69c) 



After a period of time, the Gaussian state will evolve into 
E (t) ^ U (t) E (0) (70) 
where the evolution operator 



U (t) = exp 



-i cP, 



-J—^ + V{x) + V{y)\t 
2m' / 



(71) 

Here, we assume that the detuning A is always neg- 
ative. The Schro dinger-type equation (|68|) governs the 
evolution of the wave-function of the signal light in the 
atomic medium. And the trajectory of light ray is de- 
scribed by the mean value of the coordinate operator 



Xc = (X) = / E*{t)xE{t)dx, 



(72) 




FIG. 5: (Color on /ine^Schematic illustration about the ray 
trajectory of the probe light in three dimension space. 



with X G {^1 J/j z}. An explicit calculation gives (see Ap- 
pendix B) 



Cx 



yc = ay - C 

Zc = Ct 

with the angular frequency 



1 - cos {uJqJ) 
1 — cos (uJoyt) 

m'uj^y 



i^ox = \/'^Cx2/m', LOoy = \j2(^y2/m'. 



(73a) 

(73b) 
(73c) 

(74) 



As the light travels across the atomic medium, the light 
ray - the center of the wave packet oscillates around the 
initial center in the transverse direction 



Vc 



m'oj 

J. 
m'uj 



'^Ix^ 



c 



(75a) 
(75b) 



The anisotropic motion and potential in Eq. (|68|) result 
in that, light travels in a straight line in the z-direction 
since it acts as an ultra-relativistic particle with velocity 
c, however the light oscillates in the x-y plane because it 
behaves as a non-relativistic particle with effective trans- 
verse mass m' = fc/c. If C,x2 = Cy2, the light ray is a line 
with finite length in the transverse direction. In Fig. [5l 
we schematically illustrate the wave packet center of the 
probe light in three-dimension space when Qx2 7^ Cy2- 



VII. SUMMARY 

In conclusion, we have developed a quantum approach 
for the spatial behavior of propagating light when it 
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passes through an EIT system with spatial-dependent 
external fields. By studying the dynamics of the atomic 
ensemble and the light pulse, the effective Schrodinger 
equation is derived to depict the space-time evolution 
of quasi-particles where the effective potential is in- 
duced through the steady atomic response in the exter- 
nal spatial-dependent fields. For a magnetic field with a 
spatial distribution in the transverse direction, by con- 
sidering the evolution of the Gaussian state, we showed 
that: 1) in a harmonic magnetic field, the light trajec- 
tory is a straight line; 2) in a linear magnetic field, the 
light ray bends to the direction where the magnetic gra- 
dient increases. And the deflection angle depends on four 
external parameters: the mixing angle between the sig- 
nal field and the collective atomic polarization, the wave 
number of the signal pulse, the length of the EIT gas 
cell, and the small magnetic field gradient. In an inho- 
mogeneous optical control field, we predict some novel 
results accessible on the light ray behavior. In the lin- 
ear response limit, it is found that the deflection of the 
light ray can be controlled by two controllable external 
parameters: the center position of the probe beam with 
respect to the control light, and the two photon detuning. 
In the quadric expansion of coupling amplitude, the light 
trajectory generally oscillates in the atomic medium. 

Finally we note that our study is based on a quantum 
theoretical approach. In our previous paper (l3| , the Fer- 
mat principle is applied to study the light trajectory in 
this atomic medium, and similar results are obtained, 
but that approach is a semi-classical theory , which can 
be understood in terms of the eikonal equation with the 
optical WKB approximation of our approach. Though 
this semi-classical approach can explain the most recent 
experiment about the light deflection by an EIT- 
based rubidium gas, it can not be further developed for 
the investigation of photon state quantum storage since 
the signal light was assumed as a classical field. The 
quantum approach assumes the probe light is quantized, 
thus this approach can be used to investigate the pos- 
sibility for realizing a protocol for quantum sate storage 
with spatially-distinguishable channels based on the EIT- 
enhanced light deflection. 
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APPENDIX A: FACTORIZATION OF UNITARY 
OPERATOR Uu 

Since the unitary operator Uu is exponential, factor- 
izing it means to express the exponential of a sum of 
operators in terms of a product of the exponentials of 
operators. The unitary operator Uu is an element of the 



group generated by the momentum operators Pz , Px and 
the coordinate x. Since Pz commutes with P^ and x, the 
unitary operator can first be factorized as 



where 



exp 



2m 



-Pi 



(76) 



(77) 



only contains operator P^ and x, which generate the Lie 
algebra with the basis {x,Pj:,P^, 1}. Thus, operator 
can be factorized as the form 



UL 



(78) 



and gi = gi{t) are unknown functions of time t to be 
determined. 

Mathematically, the above factorization Ansatz is 
based on the Wei-Norman algebraic theorem [l^ : if the 
Hamiltonian of a quantum system 



K 



(79) 



is a linear combination of the operators Xj that can gen- 
erate a iV-dimensional Lie algebra with the basis: 



{Xi,X2, Xk,Xk-i, ....Xn}, 



(80) 



then the evolution operator governed by H can be fac- 
torized as a product of the single parameter subgroups, 
that is , 



N 



(81) 



where the coefficients Cj(t) can be determined the "ex- 
ternal field parameters" Cj{t) through a system of non- 
linear equations. 

Now, we differentiate ((75)) with respect to t and the 
multiply the resulting expression on the right hand side 
by the inverse of ([75)) . obtaining 



2m " 
dt ^ 



^iPr. 



(82a) 



d92 „ . dgs 



dgA 
dt 



dgs 
dt 



This leads to a systems of coupled differential equa- 
tions 

.dgi 1 



dt 2m' 
■ dgs . 



dg2 
dt 

dg4 
dt 



2gii 



dgz 
dt 



dg3 



= 0, 



(83a) 
(83b) 
(83c) 

(83d) 
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The solution to these equations reads 

_ . t 
2m' 
93 = itK, 

34 = Uj. hot - 



~3 2m 



(84a) 
(84b) 
(84c) 

(84d) 



Therefore the unitary operator Uu is factorized into the 
form given in Eq. (|47p . 



II. CALCULATION OF LIGHT TRAJECTORY 
IN QUADRATIC POTENTIAL 

In the presence of the quadratic term of coordinates, 
the evolution operator U {t) in Eq. ([68]) is generated by 
the effective Hamiltonian H3 = Hz + Hy + : 



i?.^^ + y(.). 



(85) 
(86) 

(87) 



The evolution operator can be factorized as 

U (t) = {t) Uy (t) U., (t) (88) 

with 



Uz 


= Uz (t) 




(89) 


Uy 


- Uy it) 




(90) 




= {t) 




(91) 



The expectation value of the coordinator along x- 
direction is calculated as 



{x) = I E* (0) (t) xU (t) E (0) dx (92) 
E* (0) UlxU^E (0) dx 
^a,+ [ E* (0) C/t {x - a,) U^E (0) dx 



where commutation relation [x, P^] ~ iS^x is used to 
obtain the second identity in Eq. (p2|) . 

Actually, the effective Hamiltonian H^ describes a har- 
monic oscillator with its origin displaced from to other 



place. Thus, we rewrite the coordinate and momentum 
operators in terms of the creation and annihilation oper- 
ator {d^, c?}: 



V= J^ {d^ + d), 



with the inverse relation 



(dt _ d) 



d = \ I ™ 7] + ij \ Pr, 



2m'uJo 



m'uJox _ . / 1 ^ 
2 ^ V 2m' to ox 



(93) 
(94) 

(95) 
(96) 



where rj = x — and ujox is given in Eq. ([7^ . The 
effective Hamiltonian H^ can be diagonalized as 



H-r = ujnx I d^d 



by a displacement operator 
D {/3) — exp 

with 

Cxi 



f3 = p* = 



uJoxV^rn'uJo 



(97) 



(98) 



(99) 



Then the evolution operator is factorized as the product 
of three operators 

U^^D (13) e-'"^-('i"^'i+i)*D-i (/?) . (100) 

In terms of the creation and annihilation operators, the 
center of the light wave packet 



(x) =a,-H^i^ I E*{0){d^ +d)E{0)dx 



^/2m'uJox 
.sin (ujoxt) 



V2rr7uJ^ 
/3cos{u!oxt) , 2/3 



E* (0) {d^ -d)E (0) dx 



\/2m'LOox \/2m!ujo 



(101) 



Back to the coordinate and momentum operators, the 
center of the wave packet in x-direction is obtained as 
Eq. ([73k ). In a similar way, Eq. ([73b ) can be achieved. 
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